JUGGLING PROBABILITIES 



GREGORY S. WARRINGTON 



1. Introduction 

Imagine^ yourself effortlessly juggling five balls in a high, lazy pattern. Your right hand catches a ball 
and immediately throws it. One second later, it is your left hand's turn to catch and throw a ball. Then it 
is your right hand's turn again. . . . Some of your throws may be low; some high. Some balls go straight up; 
others cross over to the opposite hand. At most times, one ball lands; occasionally a hand remains empty. 
But the alternating cadence of your hands is unwavering. 

Suppose that, while you are juggling, we momentarily pause time. Certain balls are in the air — you have 
already thrown them. To avoid dropping the ball, you must make catches at certain times in the future. 
If your previous few throws were all low, perhaps you are only committed at 1,2,4,6 and 7 seconds in the 
future. On the other hand, if you had just vigorously launched a ball, you might be committed at 1,2,3,4 
and 10 seconds in the future. The set of "committed times" is your landing state ("state" for short). As you 
juggle, you wander from state to state according to what throw you have most recently made. 

Our goal in this paper is to answer the following 

Random Question. What fraction of the time is spent in any given state if every throw is chosen randomly? 

The answer, of course, depends on how we specify the following parameters: 

1. the possible states, 

2. the legal throws from a given state and 

3. the probability of making each legal throw. 

There are countless ways to make the above specifications, but we will begin our investigations with the 
model that most closely mirrors what people envision as juggling. In the next section, we describe the 
juggling universe determined by this standard model and answer our Random Question. In Section El we 
generalize the notion of a state in order to give a simple proof of our answer. Finally, in Section 01 we 
describe variations of our model. 

2. Juggling states 



Jugglers have been known to {toss, roll, spin, bomicc, drop, kick} 

their {balls, rings, clubs, flames, chairs} 

while accompanying their act with {patter, music, dance, flourishes}. 

We, however, will strip juggling down to what is (arguably) its essentials. Our juggler will juggle only 
balls and only in the air. He will not attempt jokes. Whether a hand is under a leg or behind the back 
when making a catch will not be noted. We only record the order in which the balls are thrown and caught. 
While clearly ignoring much of what makes juggling visually intriguing, our singic-mindedness will free us 
to highlight the inherently mathematical nature of juggling. 
Our juggler will: 

1. throw alternately from each hand, once a second, 

2. throw a ball instantaneously upon catching it. 



Date: February 8, 2008. 

^The author would like to thank Peter Doyle and J. Laurie Snell for their substantial help in simplifying the original proof 
of the main result of this paper. 
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3. be assumed to have always been (and be forevermore) juggling 
(this allows us to avoid starting/stopping issues), 

4. be named Magnus^. 

Every ball Magnus throws lands some (integral) number of seconds in the future. By abuse of physics, 
we will refer to this parameter as the height of the throw. Until mentioned otherwise, we now fix some 
maximum height h above which Magnus is too weak to throw. 

In the introduction, we referred to a state as a set of future times at which balls are landing. It will be 
more convenient throughout the paper to define a (landing) state v as an ft,-tuple in the alphabet {•, 
Given a state i^, we will write ut for its i-th component and write u as 



li vt — • , then Magnus will catch a ball t seconds in the future, li vt — ^, then no ball lands t seconds in 
the future. We make that convention that 1^4 = ^ for t > h. To recover our original notion of a state from i/, 
take the set {t : Vt= •}. We will denote the set of states that are /i-tuples by St}^. The subset St^ j C St^ 
will contain only those states with / ^^'s. 

As an example, consider the state f = ••^•^•••^^ G S'iio- Assuming that Magnus has just 
thrown out of his left hand, at 1 and 7 seconds in the future he will be (instantaneously) throwing/catching 
out of his right hand; at 2,4,6 and 8 seconds in future he will be (instantaneously) throwing/catching out of 
his left hand. 

When the first element of our state is v^, we will often refer to what Magnus does during the next second 
as making a "height 0" throw. Upon making a throw (even a 0), Magnus commits himself to a (usually 
different) new state u): 

If the throw is a 0, a; = ^2 ■ ■ ■ Vh ^ ■ 

where the 



1. 



2. If the throw is of height t > 1, then a; 
position (replacing vt+i)- 



V2 1^3 



Vh 



is placed in the t-th 



As an example, below we illustrate successive throws of height 2, and 6 from an initial state in which 
balls are landing at 1, 4 and 6 seconds in the future: 

2 I 1 I 1 6 



In Figure n we show all possible states and transitions when Magnus has 3 balls and is not able to throw 
anything higher than a 4. 
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Figure 1. Juggling 3 balls with maximum height 4. The edges are labeled with the throw height. 



Remark 1. In this section, we have introduced the fundamentals of "siteswap" notation — albeit in an 
untraditional form. As with many great ideas, the provenance of siteswap notation is somewhat murky. 
However, we can safely say that this notation was introduced (independently, and in various forms) in 
the mid 1980s by Paul Klimek; Colin Wright and members of the Cambridge Juggling Club; and Bengt 
Magnusson and Bruce Tiemann. For a more leisurely introduction and a bibliography of papers relating to 

After Magnus NichoUs, a juggler prominent during the early 1900s. He is reputed to be the first person to juggle five clubs. 
Reliable details of his life are outclassed by fascinating stories. 
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siteswap notation, the reader should consult 0]. General internet resources on juggling include PP and 
A recent book devoted to the mathematics of juggling is [S]. 

Both as a useful notation for juggling and for interesting enumerative combinatorics, it is preferable to 
introduce siteswap notation by defining patterns as repeating sequences of throws. Such patterns correspond 
to cycles in certain graphs; such a graph can be found in Figure ^ Many computer programs have been 
written to animate these patterns. See |3| for one cowritten by the author, or W for a comprehensive list. 

To answer our Random Question, we need an explicit description of the set of possible states Magnus 
can pass through and of the possible transitions among them. To describe these transitions, we introduce a 
family of throwing operators Oj acting on these states. Intuitively, Qq corresponds to Magnus waiting when 
no ball lands and corresponds to Magnus making a height j throw when a ball does land. Precisely, 



P2 ■ 


■ Vh ^ 


1 




if J 


= 0, 


V2 ■ 


■ • Vj+2 ■ 


■ Vh ^ 


, if J 


> 0. 



For the remainder of this and the next section, we will fix /i > and < f < h. We are now ready to 
define the directed graph which specifies the possible states for Magnus along with his possible transitions. 
Informally, the valid transitions correspond either to waiting or to making throws to empty landing times. 

Definition 2. The state graph G^j is the directed graph with 

1. Vertices indexed by St^j. 

2. An edge from u to u) whenever 

(a) a; = ©jCt') for some j with I < j < h, vi = • and Vj+i = ^; or 

(b) uj — 6o(i>') and — ^. 

The set of all edges is denoted Edges(G/i,/). 
If G Edges(G/i,/), then we will refer to i/ as a precursor of u). 

Figure El illustrates 05^2- 




Figure 2. 



Let us revisit the scenario introduced at the beginning of the paper (this time with Magnus juggling). 
When we freeze the action, Magnus is in a particular state — he is committed to making certain catches. 
If no ball is landing in one second, he is forced to wait. If a ball is landing, however, he can throw to any 
height (< h) that will not lead to two balls landing at the same time. Our assumption in the following 
is that Magnus chooses, with equal probability, one of these throws. In this manner, Magnus is taking a 
random walk on the graph Ghj. Having now phrased the question in a precise manner, we need only give 
two definitions to state the answer. 

Definition 3. The Stirling number of the second kind, |^}, counts the number of ways of partitioning an 
a-element set into b blocks (irrespective of order). 

For example, {2} = 7 as we can partition {A, B, C, D} into two parts in seven ways: 
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Table 1. Stirling numbers of the second kind {^}. 



Let V G Sth- For each integer t with 1 <t <h, define 

\0, \ivt = ^. 

Definition 4. The weight of a state i/, A(i'), is given by the formula 



(2) 



(3) ^{v) = X{{l + U^)). 

t=i 

The weight (our terminology) arises when counting permutations with restricted positions. See, e.g., [SJ 
2.4]. 

We are finally in a position to answer our Random Question. 

Theorem 5. If Magnus is juggling b balls and is able to throw them to heights at most h = b + f. then over 
the long term he spends 

A(i.) 



(4) 

of his time in the state u £ St^.f. 



1/+1J 



Example 6. Below we present the vector a whose v-th component specifies the fraction of the time Magnus 
spends in each state v of Figure |21 while juggling 3 balls at a maximum throw height of 5. 
1 
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(27,18,9,12,6,3,8,4,2,1). 



From a, we see that for not quite one third of the time, Magnus finds himself in the state corresponding to 
the usual three ball "cascade" juggling pattern. 

Remark 7. For other contexts in which the Stirling numbers of the second kind arise in conjunction with 
siteswap notation, see jHIE]- A related class of numbers, the Eulerian numbers, also arise in the enumeration 
of siteswap patterns (see II'). 



3. Proof of Theorem El 



The proof of Theorem |31 is not hard once we model the act of juggling randomly as a Markov chain: 
A discrete-time random process in which the transition probabilities depend only on the current state. To 
describe a Markov chain, we need to know the possible states f2 and the possible transitions among them. 
The latter we encode in a transition matrix P. We will denote this pair by MC(ri,P). The {i,j)-th entry of 
the matrix P is to be interpreted as the probability of transitioning from the i-th state to the j-th state in 
one step. 

In describing Magnus as taking a random walk on the states St^.f with throws chosen uniformly randomly, 
we have already described one Markov chain: MC{Sth,f, P) where P = {pv,i^)vMesth / with 



(5) 




if (i>',w) G Edges(G/i./) and vi = 
if (i>',w) e Edges(G;i,/) and vi — • , 
else. 



Suppose we have a Markov chain with states {a;^, . . . , a;''} and transition matrix P = (pij) such that the 
l-ih power P' = (pf]) of P satisfies pfj > for I ^ 0. This condition on the matrix P means that, no matter 
how long Magnus has been jugghng, he still has the opportunity to visit every state of St^j- MC{Stiij,P) 
satisfies this condition. A probability vector is a vector for which the sum of the entries is 1. The standard 
theorem of Markov chains we shall use is the following: 

Theorem 8. ( 7, 4.1.4]j Given a Markov chain as above: 

1. There exists a matrix A with lim/_^oo P' = A. 

2. Each row of A is the same vector a — (ai, . . . , ar). 

3. Each Q!i > 0. 

4. The vector a is the unique probability vector satisfying aP = a. 

5. For any probability vector tt, ttP' a as I ^ oo. 

The number ai can be interpreted as the fraction of time the Markov chain is expected to spend in the 
i-th state. Hence, Theorem |H1 tells us how to answer our question: To find the fraction of time that Magnus 
is in state w, we need only calculate the u;-th entry of the (normalized) left eigenvector for P with eigenvalue 
1! 

Example 9. We illustrate these ideas using the graph 6*4^1 given in Figure 2] The transition matrix P 



simply consists of a 1/2 for every edge except the unique edge leaving ^ • • • , which must be taken 



Notice that the sum of the entries in a row is 1. Below we list P and two of its powers. 



I ' » ' ^1 1/2 

P = I « » ^ «| 1/2 

I' ^ ' »l 1/2 

1^ • • •! 1 



P2 = I . » ^ »| 1/2 
I » ^ » »| 3/4 



1/2 

-/2 
>/4 
1/2 



I » » » ^1 17/32 

p5 = I » » ^ »| 17/32 

I » ^ » »| 17/32 

1^ • > »| 9/16 
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Already for we see that the rows appear to be converging to the same vector. And, indeed, by 
Theorem |S1 this vector is: 

a = -L(A( \. . . J ),A(r r-i— i1 ),A(rr— r-i1),A( |^. . .1 )) 

12/ 

= 1^(8,4,2,1). 
15 

Just over half of Magnus' time will be spent in the state | > » » J . In fact, this geometric distribution is 
what one would expect by examination of the graph G^^i. 

One can certainly prove Theorem [SI by guessing the vector a and then proving that it satisfies the condi- 
tions of Theorem |S| However, as is often the case in mathematics, keeping track of additional information 
enables us to write a simpler and more illuminating proof. So, we will construct an augmented Markov chain 
that is easy to analyze and from which we can obtain our original Markov chain by grouping states. 

We will augment our states by specifying which balls are landing when. To this end, define a Throwing /Landing- 



state ("TL-state") to be an /i-tuple P = Pi P2 • • • G 1,2,..., h}^ such that for each 1 < j,k,t < h, 



(6) 9t > t or 9f ^ 
and 

(7) if j ^ fc, then Vj — j ^ Vk — k. 

The first condition ensures that any ball currently slated to land has, in fact, already been thrown. The 
second condition ensures that Magnus did not have to throw two balls at once in order to get into his current 
TL-state. 

As is the case for states, we interpret = ^ to mean that no ball lands t seconds in the future. We 
interpret = i as meaning that the ball thrown i — t seconds in the past will land t seconds in the future. 
(Alternatively, the ball landing t seconds in the future will have been in the air i seconds by the time it 
lands.) A TL-state completely captures the current status of Magnus. Let Sth denote the set of all TL-states 
of length h and Stii j the subset of TL-states with / ^'s. The throwing operators Qj are extended in the 
obvious manner to act on TL-states. 

Of course, each TL-state determines some landing state simply by forgetting how long each ball has been 



in the air. So, define a map tt/i : Sth — > Sth by i^hiy) — 1^1 1^2 ■ ■ ■ Vh where 
(8) ^t^V 



if vt = ^, 

ildt e {1,2,..., h}. 



In addition to MC{Sth.f, P), we now have another natural Markov chain to consider: y[G{Sthj ,P) where 
P — iPvui) — ^-5; with 

{1, if a) = 9o(t') and vi = ^, 
jq-j-, if cD = Oj(z?) for some j with 1 < j < h, Vi = • and I'j+i = 
0, else. 

Our approach to proving Theorem [S] will be to analyze this latter Markov chain. The first step will be 
to find the vector /3 of steady-state probabilities for MG{Sthj,P). The second step will be to show that 
"lumping" states of this chain recovers our original chain MC{Sthj , P). By analyzing this lumping process 
carefully, it will be clear how to obtain the result of Theorem from our knowledge of the vector f3. 

Lemma 10. The vector of all 1 's (of length \Sthj\) is a left eigenvector for P. 

Proof. To prove this, we need simply show that the transition matrix P is doubly stochastic (i.e., that the 
rows and columns sum to 1). 

There are two cases to consider. In the first, either vt — ^ ovvt > t for every t with I < t < h. In this case, 
Magnus waited during the most recent beat. That means that the only precursor of u is 



. 1^1 • • • i^h^i 



On the other hand, il Vt — t for some t, then he just made a height t throw. Since Magnus throws at most 
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one ball at any given time, there is at most one such t. By the discussion in the proof of Lemma |14I given 
below, the previous arc of this ball could have been any of / + 1 different heights. Each of these different 
heights leads to a different precursor. Since the out-degree of any state V with Vi 7^ ^ is / + 1, each of these 
edges has weight jq-j- . This proves that the columns of P also sum to 1 . □ 

In order to fully describe the vector /3 for MC{Sth,f,P), we need to calculate \Sthj\- 
Lemma 11. \Sik.j\ = {'}Xl} ■ 

Proof. We proceed by constructing a bijection between elements of St^j and partitions with / + 1 blocks of 
the set [h + 1] := {1, 2, . . . ,h + l}. An example of the correspondence we construct is illustrated in Figure|31 

We start by associating a graph with vertex set [h + 1] to each V e Sthj. For l<t<i</i-|-l, we 
set {t, i} to be an edge of F^ if and only if there was a ball thrown h + 1 — i seconds in the past that will be 
landing t seconds in the future. In terms of 9, this means simply that i/i = h + 1 + t — i. 

Let A(P) be the set partition of [h + 1] corresponding to the connected components of Fp. Our goal is to 
show that the map V 1— > A(P) is a bijection onto partitions of [h + 1] with / + 1 blocks. 

To begin, we examine the connected components of Fj>. Suppose that {a,b}, {a,c} are distinct edges of 
Fp. As they are both incident to the vertex a, all of a, b and c must be in the same connected component. 
It cannot be true that a < b, c as then Magnus would be planning to catch two balls a seconds from now. 
Nor can it be true that 6, c < a as this would imply that a seconds ago Magnus threw two balls at once. We 
conclude that either b<a<c or c<a<b. This implies that F9 is a disjoint union of chains. 

The number of edges in Fjj is ft — / since there is an edge for each j with vj ^ ^. So, the number of 
connected components in Fjj is the number of vertices {h + 1) decreased by the number of edges {h — f). 
This yields / + 1 connected components. Hence, A(i?) is a partition of [h + 1] into / + 1 blocks as desired. 



2 3 4 5 6 7 




1 2 3 4 5 6 7 



GjTs) (46 



Figure 3. On the left we show the edges of Fp for the TL-state V = |6 ^ 46 7| . On the 
right is the corresponding partition A (9). 



Now suppose that A is a (set) partition of + 1]. We will associate a TL-state ci;(A) = cl'(A)i • ■ • w(A)/i to 
A. Consider a block {ai, q;2, am} (indexed such that tti < Qfj if i < j) of A. If rri = 1, then set ti;(A)aj = ^. 
Otherwise, set w(A)a, = h + 1 + ai — for each I with 1 < I < m. Since 1 < ai < ai+i < h + 1 for 
each /, uj{X) is certainly an ft-tuple in {^,1,2, ... ,h}. That © is satisfied follows from the inequalities 
ai < Qfi+i < h + 1. That u){X) satisfies ((JJ follows from the definition of i^(A)q, along with the fact that our 
blocks are disjoint. Hence, ti;(A) is, in fact, a TL-state. Furthermore, it can be checked that A ^ <^{X) is the 
inverse of the map u ^ A(t'). Therefore, every partition of [h + 1] into / -I- 1 blocks is obtained as a X{V) 
for a unique V G St^j. This completes the proof. □ 

Corollary 12. In the process MG{Sthj' , P), the fraction of time Magnus finds himself in state u is /h+i\ ■ 

\f+ii 

We have shown that, when juggling randomly, Magnus is as likely to find himself in one TL-state as 
another. We will now consider a new random process running in parallel with MG{Stiij, P). For v G St^j, 
define [u] — 7r^^(i/) £ Sthj. When MC{Sth.f , P) is in state 9, we define our new process to be in the state 
[7r;i(i5)]. It follows that its states are {[I'llt^esth f. This is a lumped process derived from M.G{Sthj ^P) by 
grouping together certain states. 

Certainly Magnus wanders randomly among the states of this new process. We do not yet know, however, 
that this new process is a Markov chain. In particular, it is not clear that the probability of transitioning 
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from [i/] to [a;] is independent of how long Magnus has been in [v]. If it were dependent, then the transition 
probabihtics would depend on previous states as well as the current state. Or, stated another way, the 
transition probability from [u] to [uj] would depend on which state (1 of [f] Magnus is currently in. 

For each /is [v], let Pp,,[u:] denote the probability of MG{Stfij, P,) transitioning from /i to some element 
of [u)]. Our above discussion suggests (see [3 6.3.2] for a proof) that our new process will be a Markov chain 
if the probabilities Pii.[c^] are equal for each /i G [f]. Furthermore, the transition probabilities of our 

new process will be these common values Pfi,^[u,]- 

Lemma 13. Fix f,a; € St^j- With the notation above, Pfi.[c^] = Pum for each Jl G [f]. 

Proof. If Piy.oj = 0, then certainly all of the pji^[ui] are too. Also, if Pi,^cj — 1, then Magnus must wait 
after any state /x S [f]; hence P^,[cj] = 1- So the only case to consider is when p,^^uj = JTi' ^^^^^ case, 
u) = Ojii') for some j with 1 < j < h and j^i — •. In particular, there is some j with I < j < h for which 
Vj+i = ^ and ujj = •. Suppose that Magnus transitions from /x to a particular u> £ [uj]. The previous 
equalities imply that Jij+i = <^ and ujj — j. It follows that a) is obtained from p, by making a height j 
throw. This will be one of / + 1 possible throws for Magnus as /ti ^ ^; hence he will take this choice with 
probability j^. As we have considered all cases, this completes the proof of the lemma. □ 

Define a matrix R — (?'[,y].[a;])iy,ajesth / by setting = Pi^,oj- By the above lemma, we have obtained a 

new stochastic process MC({[i/]}^g5't^ i?) that can be identified in the obvious manner with MC{Sthj,P). 
It follows from Corollarvll2l that the fraction of the time Magnus finds himself in a state v with respect to 
the process MC{Sthj, P) is TTrrrT- Hence, Theorem |5l follows from Lemma [T^ 




Figure 4. The Markov chain MC(S'i3,i, P) along with the lumped process. 

Lemma 14. 

(10) iMHA(i.)^n(i +</.,(«.)). 

t=i 

Proof. Let /? G Sthj and v G St^j such that 7r;i(j?) ~ u. Suppose vt — • and that we know the value of Vi 
for each i > t. How do these values constrain the possibilities for 9;? To start, note that the ball landing t 
seconds in the future was thrown at most h — t seconds in the past. As Magnus must have already thrown 
this ball, there are at most h — t + 1 possibilities for Vt. In addition, any ball landing more than t seconds 
in the future must have been thrown fewer than h — t seconds in the past as G Sth. Combining these 
two observations, we see that the number of possibilities for vt is reduced precisely by the number of balls 
landing after it. In particular, the possible values for Vt depend only on the Vi for i > t. 
In conclusion, the number of possible values for vt when ft = • is 

h-t + l- \{i>t:ui = •}!, 

which is just (1 + ^((i/)). Taking the product over all t yields the desired answer. □ 
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This completes the proof of Theorem |S1 

Remark 15. Consideration of the powers of the transition matrix P might lead to new asymptotic formulae 
for the Stirling numbers of the second kind: Consider Ghj and the associated transition matrix P with the 

We have already seen that the rows of P' 



last row indexed by the state uj — 



G St 



hj- 



converge to a probability vector a with a^j — Hence, for each state fi G Sthj, we get a sequence 



that converges to 



4. Other models 



By now (unless you are a very fast reader), Magnus is getting tired. He is going to drop occasionally. 
Certainly any realistic juggling model should accommodate transitions to states with fewer balls. Of course, 
if we allow Magnus to juggle indefinitely, and he drops occasionally, he will eventually run out of balls. To 
keep things interesting, we will give his assistant. Sphagnum, the option of inserting a ball into Magnus' 
pattern whenever Magnus has a wait. 

We will explore two such models in this section. In both models, Magnus can simply fail to catch a ball. 
When this happens, Magnus transitions from a state in Sthj to one in Sthj+i- The difference between 
the two models lies in what is considered to be a legal throw. In the first model, Magnus (or his assistant) 
will only throw to heights that would not lead to two balls landing at the same time. In the second model, 
however, we allow this eventuality. Of course, as Magnus is well-known to have small hands, he cannot catch 
two at once, so he will have to ignore one of the balls. Such a throw will, in terms of transitions, be treated 
no differently than a drop. 

These changes are easily incorporated into the framework we have developed. Now we define the state 
graphs corresponding to the aforementioned add-drop model (well known to college students) and annihilation 
model. 

1. Add-drop juggling generalizes standard juggling by always allowing throws of height and removing 
the restriction onvi = • for throws of positive height. (The throws of height correspond to Magnus 
dropping or waiting; those of positive height when i^i — ^ correspond to Sphagnum helping.) So, 
our graph GJ^"'' has 

(a) Vertices indexed by Sth- 

(b) An edge from u to whenever 

(i) u) = for some j , 1 < j < h such that i^j+i = ^; or 

(ii) a; = 60(1^). 
Figure illustrates G^"''. 







Figure 5. Solid arrows are regular throws; dashed are drops; dotted are insertions by Sphagnum. 
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2. In Annihilation juggling, we take the same state graph = as in the previous case. The 

difference wih be in how we assign probabihties to the edges. 
With our state graphs now defined, we can define the corresponding transition matrices. For u G St^j, 
set = \{t : 2 < t < h and i^t = ^}\- 

1. For Gl^'^, set Q = {qu,u>)u,^estu with 

1 if eEdges(Gr^), 
else. 

The denominator above comes from the fact that Magnus can either drop/wait or throw to any ^ 
occurring at V2 or later. So here, again, Magnus chooses uniformly from the available edges. 
For G^™, set R = (r.,<.).,a,GSt. with 

h-fl 





if {v,uj) e Edges(G^^ 
if (i/,a;) e Edges(G^^ 
else. 



anji ~ 
h . 
ann ' 



and 80(1^) = 1^, 
and 60 7^ i^, 



Here, you should envision Magnus (possibly with the help of Sphagnum) throwing to any height 
(< h) at random, regardless of whether or not a ball has just landed. If there is already a ball 
landing at the time he is throwing to, the latter ball gets annihilated. 

In order to describe the steady-state probabilities for the add-drop model, we must introduce the Bell 
numbers. The h-th Bell number, Bh, is defined by Bh = X]i=o {^}- ^'^^ these two models, the analogue of 
Theorem is 



Theorem 16. Let v E Sthj. 

1. In the add-drop model, Magnus spends 



B, 



h+l 



of his time in state v. 
2. In the annihilation model, Magnus spends 



h\ 



{h-f)\ 



{h + iy 



of his time in state v. 



As the denominators for our probabilities in the standard and add-drop model have natural combinatorial 
interpretations, the reader may be wondering if such an interpretation exists for the {h + 1)'' occurring in 
the annihilation model. Indeed, {h + 1)'^ counts the number of labeled rooted trees with h nodes (see JOl 
5.3.2] for details). 

Example 17. Below we give the transition matrix for Gg^'^: 



1/2 
1/2 



1/2 
1/3 
1/3 
1/2 

1/3 
1/3 



1/3 



1/3 



^. .l 








1^^^ 




1/3 








1/3 




1/3 








1/4 


1/4 


1/4 


1/4 




1/3 








1/3 




1/3 








1/4 


1/4 


1/4 


1/4 



From Theorem 1161 1. we obtain the probability vector 

1 

ct = — 
15 



(1,4,2,1,3,2,1,1) 
10 



This gives the fraction of the time Magnus spends in each of the eight possible states (the ordering of a is 
the same as that of the transition matrix) . 

As the steps in the proof of Theorem^] are analogous to those found in the proof of Theorem[Sl we do not 
detail them here. As a consolation to the reader, though, we mention two other families of juggling models 
which, although perhaps less physically realistic than the ones we have described, do lead to interesting 
answers to our Random Question. 

1. Multiplex juggling: There is no real reason (assuming you find a juggler with bigger hands) to disallow 
the catching and throwing of more than one ball at a time. The standard and add-drop models can 
both be reinterpreted with this condition relaxed. 

2. Infinite juggling: Again, given a juggler in better shape than Magnus, there is no reason to place a 
maximum on legal throw heights. In the standard juggling model, placing a probability of {p— l)/p^ 
for the throw to the j-th next available height yields tractable calculations for any p > 1. 

Rather than choosing geometric weights for the throws, one could choose probabilities from any 
convergent infinite series. For example, why not choose 6/ (tt^:/^) for the j-th available throw height? 
Or even 90/(7r'^j'*)? The sky's the hmit. 
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